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Abstract. We study the magnetisation process of the one dimensional spin-1/2 antiferromagnetic Heisen- 
berg model with modulated couplings over j = 1, 2, 3 sites. It turns out that the evolution of magnetisation 
plateaus depends on j and on the wave number q of the modulation according to the rule of Oshikawa 
et al. A mapping of two- and three- leg zig-zag ladders on one dimensional systems with modulated cou- 
plings yields predictions for the occurence of magnetization plateaus. The latter are tested by numerical 
computations with the DMRG algorithm. 
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1 Introduction 

The formation of gaps and plateaus in the magnetisa- 
tion process of one-dimensional (ID) spin-1/2 antiferro- 
magnetic Heisenberg models has been studied intensively 
during the last years |,|,il|,i||,|,UlO|. 

Oshikawa, Yamanaka and Affleck |l|pointed out the 
crucial role which play the soft modes, predicted by the 
Lieb-Schultz-Mattis (LSM) construction . For example 
in the case of the ID spin-1/2 Hamiltonian with nearest 
neighbour couplings and a homogeneous external field B: 



H(B) = Hi-BS3(0), 

N 



2^ Si • Si+j , 
(=1 



1,2,.. 



N 

S,(<z) = ^e*''5r, a = 1,2, 3, 
1=1 



(1.1) 
(1.2) 

(1.3) 



the ground state |0) has momentum ps = 0, tt and total 
spin Si^ = St ~ NM{B), where M{B) is the magnetisa- 
tion. 

The LSM construction (ll],|l^ leads to gap-less excited 
states \k): 




= U'-'IO), U = exp 

with momenta 

Pk=Ps + kq3{M), q3{M) = Tr{l-2M), 



(1.4) 



(1.5) 



e.g. for M — 1/4 one finds a four fold degeneracy of 
the ground state with momenta pk — Ps + fc7r/2, k ~ 



0, 1, 2, 3. The magnetisation curve for the system ( |l.l| ) has 
no plateaus. The latter app ear, if translational invariance 
is broken, by adding to (1.1) a periodic perturbation with 
wave vector q, e.g.: 



I)j{q)^^[r>j{q)+-D,i-q)], (1.6) 

N 

B,{q)^2j2e"'%-Si+j^ J = l,2,.... (1.7) 



So far the case j — 1 has been studied in detail ||l^. A 
pronounced plateau has been found in the magnetisation 
curve if the period (^/tt = 1/2,1/6,1/3) coincides with 
the first soft mo_de {k ^ 1), q = q3{M) = 7r(l - 2M), i.e. 
a perturbation Di(7 r/2) generates a plateau at M = 1/4. 
However, adding to ( |l.l| ) a perturbation Di(7r) no plateau 
has been observed at M = 1/4. In the latter case the sec- 
ond soft mode at 9 = tt (fc = 2) coincides with the period 
of the perturbation. Therefore, the question arises, why 
certain possibilities for the formation of plateaus, which 
are allowed according to the quantisation rule of Oshikawa 
et al. 0], are nevertheless not realized with a given per- 
turbation. One possible answer to this question has been 
given in Ref. [|l^: The efficiency of the mechanism to 
generate a plateau, by means of a periodic perturbation 
Di(g), crucially depends on the magnitude of transition 
matrix elements {n\Di{q)\0) from the ground state |0) to 
the low-lying excited states |n). For example for M = 1/4 
the transition matrix elements (n|Di(7r/2)|0) turn out to 
be large whereas (n|Di(7r)|0) are small. The effect can 
be seen directly in the static dimer-dimer structure factor 
(0|Di(g)Di(-g)|0) for M = 1/4 (cf. Fig. 4 in Ref. [0), 
which has a pronounced peak at q ~ tt /2 but no peak for 
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g = TT. Of course all these stat eme nts only hold for the 
nearest neighbour Hamiltonian (1.1). 

Indeed, Totsuka [|| has recently observed a magnetisa- 
tion plateau at M — 1/4, which was created by adding to 
the Hamiltonian (|l . l| ) a perturbation Di(7r) and a strong 
next-to-nearest neighbour coupling. 

In this paper we show that the magnetisation plateau 
at M = 1/4 can be realized as well with perturbations 
of period q = tt, if the perturbation operator is chosen 
properly. We will discuss i n d etail the situation with the 
operator D2(g) defined in (|l.(j). Note, that both operators 
Di(g) and T>2{q) have the same momentum- and spin- 
symmetry properties; they change the momentum by ±q 
and do not change the total spin. They only differ in the 
isotropic spin-spin couplings, which extend over nearest 
neighbours in T)i{q) and over next-nearest-neighbours in 

D2(g). 

In Sec. ^ we compare the static structure factors, Sj ^ 
(0|Dj(g)Dj(— q)|0), j — 1, 2, in the presence of a magnetic 
field with magnetisation AI — 1/4. In Sec. ||we study the 
magnetisation plateaus induced by the periodic perturba- 
tion D2('z) at q = TT. 

The zig zag ladder with two (and three) legs - recently 
investigated in Ref. Q - can be mapped on a ID sys- 
tem with translation invariant coupling over one and two 
(three) lattice spacings and a translation invariance break- 
ing coupling of the type I)2{q)- In Sec. ^ we analyse the 
sequence of magnetisation plateaus, which appears in the 
zig-zag ladder system. 

It should finally be added that the operators and Hamil- 
tonians refer to periodic boundary-conditions in leg direc- 
tions. The DMRG results given in Sec. ^, however, have 
been obtained using open boundary-conditions along the 
legs. 



2 Signals of soft modes in static structure 
factors 



In this section we discuss some properties of the static 
structure factors {0\'Dj{q)'Dj{—q)\0), j = 1,2, of the Hamil- 



The operators Dj(g), j — 1,2, defined in Eq. (1.7) 



tonian (1.1). In order to compute static structure factors 
we use periodic boundary conditions and exact Lanczos 
diagonalizations up to iV = 24 sites. 

Soft modes, as predicted by the LSM construction, can 
be seen directly as zeros in the dispersion curve [0 : 

uiq, M) = E{ps +q,S+l)- E{ps,S), M = SI./N, 

(2.1) 

where E{p, S) are the lowest energy eigenvalues with total 
spin S = St and momentum p. The ground-state momen- 
tum ps in the sector with total spin S is known to be or 
7r|ll^. The zeros of (2J) - in the limit N ^ oo - appear 
at the soft mode momenta q = q^''\M). For example, at 
M = 1/4 three zeros at q/n — 0, 1/ 2, 1 emerge in the 
dispersion curve for the Hamiltonian ( [l.l[ ) (cf. Fig. 3 of 
Ref. M). 



and the dis- 



commute with the total spin squared 
persion curve (2.1) describes the lowest-lying excitations, 
which can be reached with these operators. The transition 
matrix element 



(n|D,(g)|0), j = l,2, 



(2.2) 



from the ground state |0) with total spin S = St and 
momentum p^ to the excited states \n) with momentum 
p ~ Ps + q enter in the corresponding static structure 
factor: 



5,(g,M)^l^(l-MI(HD,((?)|0)| 



(2.3) 



The signals of the soft modes in the static structure fac- 
tor therefore measure the magnitude of the transition ma- 



trix elements (2^). In Figs. yT|(a) and 2T(b) we compare 
the g-dependence of the static structure factors Sj{q, M = 
1/4). Indeed we observe remarkable differences. We find a 




Fig. 2.1. C omp arison of the q-dependence of the static struc- 
ture factors (2.3) at M = 1/4 for the dimer operators (1.6) over 
nearest [ (a) j = 1] and next-nearest - nei ghbour [ (b) j — 2] cou- 
plings in a model with Hamiltonian (1.1 ). The different symbols 
belong to system sizes: TV = 20 (■), 16(o), 12 (•). 



pronounced peak in both structure factors Sj{q, 1/4), j — 
1, 2 at the soft mode q^^-'(l/4) = 7r/2. The size of the peaks, 
however, differ by an order of magnitude. There is no peak 
at the second soft mode g'^^^(l/4) tt. Note the different 
behaviour of the structure factor. For q ^ 0: 5*1 (g, 1/4) 
converges to zero whereas £'2(9, 1/4) approaches a max- 
imum < 0.3 . This feature will play a crucial role, if we 
add to (IT) a periodic perturbation S ■ Dj(7r), j = 1,2 



of strength 6. The perturbation is invariant under trans- 
lations by two lattice spacings. Eigenstates of the per- 
turbed Hamiltonian are constructed by a superposition of 
momentum eigenstates with p — and p = tt. The re- 
duction of the BriouUin zone is taken into account in a 
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modification of tlie dimer operators: 

N/2-l 



Tlie corresponding structure factors, defined by: 




(2.4) 



S^(q,M) ^ 1 [(0|D^-(cz)D^-(-q)|0) - |(0|D^«(g)|0)| 

(2.5) 

are symmetric under the mapping q ^ n — q. 

In Fig. |]^(a) we present the static structure factor 
Si{q, 1/4), K = 0, 1 - obtained from the ground state of the 
Hamiltonian (1.1) with perturbation 6-'Di{tt), S = 0.5. In 
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Fig. 2.2. Comparison of the g-dependence of the static struc- 
ture factor ( p.5| ) at Af = 1 74 f or the dimer operators (2.4) in a 
model with Hamiltonian (1.1) and perturbation 5 ■ (tt), k = 
0(»), 1(0); j = 1(a), 2(b). The perturbation strength is 5 = 0.5, 
the system sizes are = 20, 24. 



We dicsuss the same Hamiltonian for which we have com- 
puted the static structure factor in Sec. 2.2. For this reason 
we also apply periodic boundary conditions in this section. 

The possible position of the plateaus is given by the 
quantisation rule of Oshikawa tt al. However, whether 
or not a plateau really appears, crucially depends on the 
type of the pert urbation operator. As an example we com- 
pare in Fig. the evolution of the magnetisation plateaus 
generated by the operators Di(7r) and D2(7r). 

The magnetisation curves - generated with Di(7r), cf. 
Fig. 



3.1 



column (a) ~ show a plateau at M = 0, rapidly 
increasing with the strength d of the perturbation. This is 
the well known gap induced by dimerisation. There is no 
plateau at M = 1/4 in the whole 6 range {0 < S < 0.7). 

In contrast the magnetisation curves - generated with 
D2(7r), [cf. Fig. 3.1 column (b)] - have no plateau at M = 
0. For S > 0.4 a plateau appears at M = 1/4. A finite-size 
analysi s of the plateau width yields the (5-evolution shown 
in Fig. 



3.2 



Note in particular, the drastic change in the plateau 
width at (5 = 0.7. This feature is associated with a change 
in the ground-state quantum numbers in the fixed S^- 
sectors. For S < 1/2 all ground states (0 < M < 1/2) have 
the standard momentum 0, tt. The first change happens in 
the sector M — 1/2—1 /N at (5 = 1/2, where the ground state 
is degenerate with momentum p = 0, tt and p = ±7r/2. 

For larger values oi 6 {6 > 1/2) we observe a different 
ground-state behaviour in the sectors with: 

0<M<MoiS), and Mo((5) < Af < 1/2. (3.1) 

In the first regime of ( ^j| ) the ground state has still mo- 
mentum p — 0, TT. In contrast, in the second regime of (3.1) 
the ground-state momentum alternates between p = 0, tt 
and p — ±7r/2. The magnetisation Mq{S), which sepa- 
rates the two regimes passes the plateau M = 1/4 exactly 
at 5 = 0.7. 

The plateau structures in the magnetization curves of 
chains with periodic perturbations can be seen already on 
rather small systems {N w 24). 



both structure factors we only find a peak at g = 7r/2 but 
no peak at g = 0, tt. The situation for the static structure 
factors S^{q, 1/4), k = 0, 1 is shown in Fig. U(b). Here 



the ground state has been computed for Hamiltonian (1.1) 
with a perturbation S ■ D2(7r), 6 = 0.5. Note the different 
behaviour of the two structure factors 52(9,1/4), k = 
0,1. 52(9,1/4) has its maximum at q = 0, tt, whereas 
5|((7, 1/4) has a maximum at q = tt/2. We therefore ex- 
pect that the operator D2(7r) generates a plateau at Af — 
1/4, whereas the operator Di(7r) does not. 



3 Magnetisation plateaus induced by periodic 
perturbations 

Periodic perturbations of the type D, (g) - added to the 
nearest neighbour Hamiltonian ( |l . l[ ) - generate a charac- 
teristic sequence of plateaus in the magnetisation curve. 



4 Magnetisation plateaus in zig-zag ladders 

The magnetic properties of zig-zag ladders have been re- 
cently investigated by Cabra, Honecker and Pujol [|3) by 
means of bosonization techniques and numerical analysis. 
In order to understand the appearance of plateaus in the 
magnetisation curve, a mapping of the ladder system onto 
a ID spin chain with couplings over short range distances 
is useful. For normal spin ladders with / legs this mapping 
leads to couplings over nearest neighbour and to couplings 
over / lattice sites Q . The couplings over I lattice sites ap- 
pear along the legs, whereas nearest neighbour couplings 
form the rungs of the ladder. At the endpoints of the rungs 
nearest neighbour bonds have to vanish to avoid the ap- 
pearance of diagonal couplings in the ladder. A Fourier 
analysis of the translation invariance breaking terms: 



(4.1) 
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10 0.5 1 1.5 



Fig. 3.1. Comparison of the plateau evolution in the magneti- 
sation curves of a model with Hamiltonian (1.1) and perturba- 
tion 5-t)j{Tr),j = 1(a), 2(6), S — 0.5, 0.6 (from top to bottom). 
The solid lines represent the midpoint magnetisation curves, 
together with extrapolated values of the upper and lower crit- 
ical field B'^ and B^, respectively, deduced from system sizes 
iV = 8, 12, ... ,24. 




In the following subsections we consider different re- 
alizations of zig-zag spin ladders, which have been dis- 
cussed recently by several authors |l^,|l5|,|l6). We map 
those systems to ID systems with appropriate couplings 
over j = 1, 2, 3 neighbours. For some of the following sys- 
tems we expect more than one plateau in the magnetiza- 
tion curve, therfore we use in this section DMRG calcu- 
lations, which also means we apply open boundary con- 
ditions, to get a finer resolution, i.e. more steps in the 
finite-size magnetization curves. 



4.1 The two-leg zig-zag ladder 

The mapping of the two leg zig-zag ladder onto a ID sys- 
tem with short range couplings is shown in Fig. 4.1. The 



_2 J2(i-S) _4 j2(i+s) 




Fig. 4.1. Mapping of the tw o-leg zig-zag ladder onto a ID 
system with Hamiltonian (4.4). The solid lines represent the 
nearest neighbour interactions with coupling Ji and the dashed 
lines the next-nearest-neighbour interactions with couplings 

J2(l±5). 



Hamiltonian can be written as: 

H = ,/iHi + J2(l + S)n+s + J2(l - 5)H_5, (4.4) 

N 

Al ■ S4;±2 + 84; + ! • S4;±2+l] • (4.5) 



1=1 



The Fourier analysis (4.1) of the translation invariance 
breaking terms yields in this case: 



Fig. 3.2. The (5-evolution of the M = 1/4 plateau width 
A{5) — B^ — B^ in a model with Hamiltonian (IT) and 
a perturbation 5 ■ 02(71), deduced from system sizes A*' = 
12,16,20,24. 



leads to a prediction of magnetisation plateaus at: 



7 1 Z 

M ^ MF = , 



(4.2) 



where Z is integer and runs over the sequence 



I 1,2,3,... ,Z/2 :/even 
[1,2,3,... ,(?-l)/2 :/odd. ^ ■ ' 

The prediction is based on the assumption that a plateau 
only appears if one of the wave vectors q in the Fourier 
analysis (4.1) of the perturbation coincides with the first 
soft mode [k = 1). 



H 



i=i,2 



N 



Z_,COS - ^ ) • S/+2. 



(4.6) 



Therefore we expect magnetisation plateaus, if the wave 
vector q/tt = 1/2,3/2 meets the first and second soft 
mode: 



q = 



2 
37r 



Ml 



Ml 



1 1 

2 ^ 4fc' 

1 3 

2 " 4fc' 



forfc = l,2, (4.7) 
for fc = 2. (4.8) 



We have looked in particular for plateaus at Af| = 1/8 
and Mg =3/8 induced by the second soft mode fc = 2. 
The situation at J i = 1, J2 = 2 and 5 = 0.6,0.8,1.0 is 
shown in Fig. |4.2| , where the emergence of a plateau at 
AI ^ 1/8 is visible. The effect disappears if we change the 
ratio a — J1/J2 in both directions a < 1/2 and a > 1/2. 
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Fig. 4.2. Magnetisation curve of the two leg zig-zag ladder 
with Hamiltonian (4.f) with couplings: Ji — 1, J2 = 2, S — 
0.6, 0.8, 1.0. The system size is iV = 48. 



4.2 Three-leg zig-zag ladder with periodic rung 
couplings 

The mapping of the three-leg zig-zag ladder onto a ID 



system with short range couplings is shown in Fig. 4.3 
The corresponding Hamiltonian can be rewritten as: 




Fig. 4.3. Mapping of the three leg zig-zag ladder with peri- 
odic rung coupUngs onto a ID system with Hamiltonian ( L^ j . 
The solid lines represent nearest neighbour interactions with 
couplings Ji, the dashed lines next-nearest-neighbour interac- 
tions with couplings J2, and the dashed-dotted lines interac- 
tions over third neighbour interactions with couplings J3. The 
dotted lines denote interactions with nearest neighbour cou- 
plings Ji , as they are enforced by periodic boundary conditions 
along the rungs. 



We have computed the magnetisation curves with pe- 
riodic rung couplings and open boundary-conditions along 
the legs for the following values: 



3 3 



(4.11) 



and system sizes: N = 3 x L, L = 8,12, 16, 20, 24. It 
turns out, that the midpoint extrapolation a la Bonner 
and Fisher [0 still shows a surprisingly large finite-size 
effects for the ope n bo undary conditions. However, as can 
be seen from Fig. iA, there is still a clean signal for two 
plateaus at M — 1/6 and M = 1/3. 

In Ref. p^ , the magnetisation c urve has been com- 
puted for the same set of couplings (4.11), but with dif- 
ferent boundary-conditions along the rungs, which they 
call periodic boundary-conditions of type A, B or C. No 
plateau at all is found for type B, one plateau at M =1/6 
is found for type A and C. This is a first indication that the 
formation of plateaus critically depends on the boundary- 
conditions along the rungs. 



CQ 




B 

Fig. 4.4. Magnetisation curve of the three leg zig-zag ladder 
with periodic rung couplings and Hamiltonian (j^^ with cou- 
plings: Ji — 1.5, J2 = 1.5. The system size is A'^ = 72. The 
solid line is the midpoint magnetisation curve. The behaviour 
at the end of the plateaus is extrapolated from data of sys- 
tems A'' = 24, 36, . . . , 96. The finite-size effects are largest in 
the region around the upper critical field of the plateaus. 



H 



N 

1=1 



27r 



n s, • s 



1+2 



(4.9) 



The last term on the right-hand side breaks the transla- 
tion invariance of the ID system and we therefore expect 
magnetisation plateaus for 



1 



1 



Mi ^ -, Ml = -, 
3 g, 63. 



(4.10) 



if the wave vector q = 27r/3 meets the first (fc ~ 1), second 
(fc — 2) and third (fc — 3) soft mode, respectively. 



4.3 Three-leg ladder with open rung couplings 

The mapping of this system onto a ID syste m w ith short 
range couplings can be seen again from Fig. 4.3. We only 
have to remove the dotted nearest neighbour bonds, which 
implement the periodic rung couplings. This changes the 
Hamiltonian in the following manner: 



+ J3H3. 



N 



Hj — 2 cos 



1=1 



+ 2-j) 



S/-S 



1+3 



(4.12) 
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The two Hamiltonians (4.9) and (4.12) differ in tlie relative 
weiglit of the translation symmetry breaking terms. 

It was found already in Ref. |Q that for the couplings 

([4.1l|) (Ji = 3/2, J2 = 3/ 2, J3 = 1) the magnetisation 



curye for the Hamiltonian (4.12) evolves two pronounced 
plateaus at M = 1/6 and at M ~ 1/3 generated by the 
first and second soft mode, respectively. This is confirmed 
in our computation. Switching from antiferromagnetic to 
ferromagnetic leg couphng ( J3 — —1) we find spontaneous 
magnetisation at M =1 / 6 whereas the plateau at M = 
1/3 disappears (Fig. [4. 51). 



CO 




Fig. 4.5. Magnetisation curve of the three leg zig- zag l adder 
with open rung and leg couplings and Hamiltonian ( [4.12| ) with 
couplings: Ji = J2 = 3/2, J3 = —1. The solid line represents 
the midpoint magnetisation, extrapolated from system sizes 
A'' — 24 — 72. The thin line is the magnetisation curve for 
N = 48. 



4.4 The Kagome like three spin ladders 

This system has been studied recently in Ref. ||l^ for rea- 
sons which we expl ain b elow. Its couplings are defined by 
removing from Fig. 4.3 the middle leg and the dotted ver- 
tical lines. The mapping onto the ID system with short 
range couplings then leads to the Hamiltonian: 



H — - ^ J^Hr 



r- 




4 








3 






N 


J2 


E 




1=1 




N 



N 

E 



J3 ^ cos 

1=1 



cos 



3 ^ 



S, ■ S 



i+i- 



I' 



27r 



;+3, 



(4.13) 



Equation ( 4.13 ) differs from Eq. ( 4.12 ) in the additional 
translation symmetry breaking term over three lattice spac- 



ings. The Kagome like three leg ladder is interesting be- 
cause it shares with the two dimensional Kagome lattice 
the property that there is a high density of low-lying sin- 
glets. A different Kagome like three leg ladder system has 
been analysed recently by Azaria et al. . 

The Kagome lattice itself seems to have a singlet-triplet 
gap. The authors of Ref. ||l5| tried to find out whether 
such a gap (i.e. a plateau at M = 0) exists as well for the 
Kagome like three spin ladder. Extrapolation from DMRG 
results for systems up to 120 sites show that the system 
is gap-less for values of the leg spin coupling J in the 
interval 0.5 < 2J3 < 1.25, Ji = 1/2. Note, that the trans- 
lation invariance breaking terms in (4.13) do not generate 
a plateau at M = but at M = 1/6 and AI = 1/3. 

Our results for fixed couplings Ji — J2 — 1-0 and in- 



creasing couplings J3 — 0.4,0.6, 1.0 are shown in Fig. 4.6 
For J3 = 0.4 we find spontaneous magnetisation at M = 
1/6. This phenomenon is a consequence of the Lieb-Mattis 
theorem , as was pointed out by the authors of Ref. ||l^ . 
In condensed matter physics this phenomenon is called 
ferrimagnetism. The magnetisation curve for J3 = 0.6 
starts with a steep increase at _B = and reaches quickly 
the plateau at M = 1/6 for a small value of i? « 0.1. Going 
to larger values of the coupling J3 — 1.0, the slope of the 
magnetisation curve and the plateau width at AI = 1/6 
are reduced. At the same time, we expect the appear- 
ance of a small plateau at M = 1/3. Beyond this point 
the system jumps into the saturation magnetisation. This 
phenomenon is called meta-magnetism E0| . 




B 



B 



B 



Fig. 4.6. Magnetisation c urve of the Kagome like three spin 
ladder with Hamiltonian ( 4.15 ) and couplings: Ji = J2 ~ 
1, J3 = 0.4,0.6, 1.0. The system size is iV = 48 for (a) and 
(b) and due to the larger finite-size effects N = 72 for (c). The 
dashed lines show the extrapolted postions of the plateau ends 
deduced from Lanczos diagonalizations with periodic boundary 
conditions. 



We finally want to give some comments on the strong 
finite-size behavior of the presented DMRG calculations 
of magnetization curves especially appearing at the up- 
per critical fields of the shown magnetization plateaus. 
Reanalyzing the 3-leg zig-zag and Kagome like ladders 
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of this section for system sizes N — 12, 18, 24 and peri- 
odic leg boundary conditions applying standard Lanczos 
techniques led to a straight affirmation of the presented 
extrapolations. It moreover showed that the strong finite- 
size effects at the upper plateau edges are no genuine ad- 
ditional features but a simple consequence of the open leg 
boundary conditions used for the DMRG calculati ons (see 
e.g. plateau boundaries additionally given in Fig. 4.6). In 
addition, the evaluations with periodic leg boundary con- 
ditions could make clear the existence of an additional 
plateau (m = 1/3) for the Kagome like 3-leg ladder shown 
in Fig. 16 (c), i.e. J3 — 1.0. Again, at the upper critical 
field strong finite-size corrections hinder the identification 
of the m = 1/3 plateau on the basis of the shown DMRG 
results. For the two other cases -J3 = 0.4, 0.6- additional 
plateaus could be excluded. It remains to be summarized 
that evaluations with both types of leg boundary condi- 
tions agree -if possible to obtain- in the thermodynamic 
limits of the considered quantities. Periodic boundary con- 
ditions, however, show much smaller and more controlled 
behavior of finite-size effects while the open boundary 
conditions require considerably larger system sizes for an 
equal quality of extrapolation. 



once. The links which do not lie on the ring, define further 
reaching and modulated couplings. However, these map- 
pings become unique, if we postulate that the range of 
the couplings is minimal. This means for normal ladders 
with /-legs, that the corresponding ID system only con- 
tains modulated nearest neighbour couphngs [Di((7)] and 
(translation invariant) couplings over I lattice sites. Here 
magnetisation plateaus appear if the wave vectors q of 
the magnetisation plateaus meet the first soft mode. The 
situation is different if we map zig-zag ladders on quasi 
ID systems. Modulations of the next-nearest-neighbour 
couplings [D2(g)] emerge, which generate magnetisation 
plateaus if the wave vector q of the modulation meets ei- 
ther the first or the second soft mode. 

Finally, we studied the magnetisation process of three 
leg zig-zag ladders with various boundary-conditions along 
the rungs. The boundary-conditions change the weight of 
the terms which modulate the couplings over one, two 
and three lattice spacings in the ID Hamiltonian. This 
again affects the formation of plateaus at M = 1/6 and 
AI ~ 1/3, respectively. 
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5 Discussion and Conclusions 

The Li eb-S chultz-Mattis construction of gap-less excited 
states (1.4) in quasi one-dimensional spin- 1/2 quantum 



spin systems demands translation invariance and short 
range couplings. To our knowledge systems which satisfy 
these conditions have no plateaus in their magnetisation 
curve for Af > 0. In this paper we have studied the effect 
of a modulation of the couplings over j — 1,2,3 sites [cf. 
Eq. (1.6)] on the magnetisation process. 

According to the quantisation rule of Oshikawa, Ya- 
manaka and Affleck, we expect plateaus if the wave num- 
ber q of the modulation coincides with one of the momenta 
of the soft modes g^'^) = fc7r(l - 2M), fc = 1, 2, . . . . We 
have found, that the modulation of the nearest neighbour 
coupling [with the operator Di((7)] generates one plateau 
at M = 1/2 — q/iT, i.e. if q coincides with the first soft 
mode (j(i)(Af). 

For q = TT, the modulation of the nearest neighbour 
coupling with the operator Di(7r) leads to the well-known 
singlet-triplet gap, i.e. a plateau at Af = - which opens 
with the strength 6 of the perturbation as S'^^^ • 

The modulation of the next-nearest-neighbour coupling 
with D2(7r) does not affect substantially the magnetisation 
process if (5 < 0.4. There is neither a plateau at Af = 
nor at M — 1/4. For S > 0.4, however, a plateau opens 
rapidly at A/ = 1/4 and shrinks again for (5 > 0.7. We 
have found that this effect is correlated with a change in 
the quantum numbers of the ground state. 

Spin ladders with I legs can be mapped on one dimen- 
sional systems with modulated short range couplings. Of 
course these mappings are not unique. There are many 
possibilities to put a ring with nearest neighbour cou- 
plings on a ladder in such a way, that each site is passed 
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